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Abstract 

Ragistration of tharaal iaagas is coaplicatad by dlstinctlva diffarancas 
in tha appearanca of day and night features needed as control in the registra- 
tion process. These changes are unlike those that occur between Landsat 
scenes and pose unique constraints. Experiaentation with several potentially 
proaising techniques has lead to selection of a fairly siaple scheme for 
registration of data froa the experiaental thermal satellite HCMM using an 
affine transformation. T%io registration examples are provided. 


Introduction 

The technique of theroial-inertia mapping (Watson, 1971; Watson and oth- 
ers, 1972; Pohn and others, 1974; Watson, 1975; Kahle and others, 1976; Price, 
1977; Pratt and Ellyett, 1979; Watson, 1981b) requires co-registerlng of 
images acquired at different tines to construct temperature-difference and 
theraal-lnertia images. In addition, refinements to use topographic correc- 
tions (Cillespie and Kahle, 1977; Hiller and Watson, 1977; Watson, 1981a) and 
to coapare theraal data with other data bases, such as Landsat images and 
geophysical and geologic naps, require further registration to a topographic 
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b«t«* Thcnal iMg«s havt mm imtqM eharMCcrittlet that aak* this rtgls- 
tratiea procsss oors esoptn than rstiscration of stsnisH .phototraphs and 
Landaat ioagot. Tha Inharant tMporal bshavior of tha atirfaca taaparature, 
‘that dapandt <m albado, tharaal and aoistivc pro^rtias, and slops, alavation, 
and oolstara affects, aokas tha scene appear substantially different on day 
and night scenes. Texture and topography often are portrayed in an entirely 
different manner on theraal images acquired at different times during the day 
and night (Sabins, 1969; Rotmn and others, 1970; Offield, 1975). The selec- 
tion of reproducible control points associated with various features thus can 
be a difficult task aiul one subject to substantial error. For example, mislo- 
cation of a feature on an image with respect to a corresponding topographic 
ridge can result in a topographic correction being applied with the wrong 
sign! The Heat Capacity Mapping Mission (HCMM) data also presents an addi- 
tional problem because the 500-m digital resolution masks many distinctive 
features— primarily cultural— that can be seen on the 80-m resolution Landsat 
data. 


Errors resulting from the misregistration between images and among images 
and various data bases have been discussed previously (Hiller and Watson, 
1977). The purpose of this brief paper is to describe our current experie.ice 
with the registration of HQOt data and to suggest the use of a simple method 
for registering the data that we have found to be more accurate than the 
registered products currently being made available. 
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Method 


Ue experimented with Iterative registration^ with uping drainage- 
intersection asps for control, and with cross correlation techni<)ues but none 
were found to be satisfactory. The procedure finally used to register the 
image pairs was to select control points and to map the night thermal image to 
the daytime thermal and reflectance images using an affine transformation on a 
1300 by 1100 pixel image. The resulting image re. 'stration was accurate to 
better than two pixels (rms) and does not exhibit the significant misregistra- 
tion that we have noted in the temperature-difference and thermal-inertia 
products supplied by NASA. 

The affine transformation was determined using simple matrix arithmetic— 
a step that can be performed rapidly with matrix hardware on a minicomputer 
such as a Hewlett-Packard HP98A5^^. 

Consider a set of control points in the day image and the corre- 
sponding points Xj*,yj^’ in the night image. A best-fit affine transformation 
can then be computed: 

let S' ■ (xj^.yj^) 

S - (x^.y^.l) 


— ^Use of trade names is for descriptive purposes only and does 
not constitute endorsement by the U.S. Geological Survey. 
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then the desired affine traneforaation 



tatisfies the oatrix equation S* ■ ST and thus 

T ■ (S^S)”^ (S^S*) vhere is the transpose of S. 

Results 

The affine transformation provides a rotation correction for the Inclined 
day and night orbital tracks (see figure 1), an origin shift, and a scale 
change. Our previous experience with aircraft data suggested the need to 
correct for platform instabilities due to pitch, yaw, and roll of the aircraft 
and this has required non-linear or piece-vise linear corrections. Ue did not 
find this necessary with the two HCMN image data sets that we have analyzed. 

Presumably the satellite, being more stable, is l*ss subject to these effects. 

The two most important considerations in selecting control features are 
that they are readily identifiable on both day and night data and that their 
locations cover the scene. We experimented with selecting different types of 
features as control points. The best features— those most reproducible on 
both day and night data— were the water-dam interfaces of reservoirs. Dis- 
tinctive outlines of other water bodies were also reasonably identifiable. 
Drainages and topographic features were found to be the least reliable due to 
changes in their day-night appearance. In the first registration example 
(figure 2), nine water body features were used as control. In the second 
example (figure 3), few water bodies exist, and only three of the six selected 
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control features were water-aaaoclated. The other three pointe are topograph- 
ic features. The accuracy of identifying these features was in part estimated 
by the magnitude of the residual errors. A high residual was indicative of 
Ahe feature being misldentified. Figure 4 shows an example of the control 
vectors (a) and of the residuals (b) from the Powder River Basin scene regis- 
tration. 

Ve discovered that a simple transformation was satisfactory over large 
parte of an entire image in two relatively arid regions of the western United 
States, using a few accurately determined control points at bodies of water. 
The procedure should be more easily applied in other areas where bodies of 
water are likely to be more abundant. A listing of the computer program to 
perform this registration is provided in Appendix 1. 

Automatic registration using cross-correlation failed when we attempted 
to apply it to one site because of a pronounced topographic grain that pro- 
duced a significant error along the direction of the grain. 

Some sample results of our registration efforts are provided to illus- 
trate the method and its simplicity. Figure 2 shows the transformed night 
image compared with the day image and the original untransformed night 
image. The second example is the Cabeza Prieta, Arizona, area in figure 3. 
Day, night, and thermal-inertia Images are presented, the latter being comput- 
ed from the registered day thermal, night thermal, and reflectance images 
using a new algorithm (Watson, 1981b). Figure 4 presents a comparison between 
the original control vectors and their residuals after applying the method to 
HCMM data of the Powder River Basin, Wyoming, area. 
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HtooRAn cconx4 

c..... 

c 

c mnoTE ccMtiNO mruvf pkocedurcs 

C..... U. «. OCOCOOICAL £l«iVEV. CCNVCA. COLOAACIO 

C BRANCH OP PETAOH4Vfr2C£ AND P£nOTE SENSINO 

• C DON L. SAUAT2KV 

C* • • • • 

c 

C 0E0HX4 OEICAATES A RECTIFIED IMAC€ FILE FAOH A DISTORTED INAOE 
C FILE ANIi FROM COEFFICIENTS FCM RECTlFICATim. R AND S. DETERHINED FOR AN 
C AFFIIC TRANSFORTMTION OF THE DISTORTED FILE. 

C INPUT FILE STRUCTU^tE CC<NSISTE CF A »«AD£R RECORD CONTAINING TUO 
'C INTEGERS FC«( LINE LENGTH. LENREC. IN PIXELS AND NUMRER OF Llf^S. 

C NC4TEC6. NW ECS NUNRER CF DATA RECCM%DS FOLLOW. EACH RECORD CONTAINING LENREC 
C BYTES OF 8-1 IT DATA. C*UTF-UT FILE »«ADER RECORD CONTAINS TWO INTEGERS tF LINE 
C LENGTH. NPXOUT* ANt< NIA«ER OF RECi^DS. U1AX. INPUT PARAHETERS. IPIXIN AND NPI 
C XIN. ALLOW TAKING A SUSSET OF THE INPUT FILE. OUTPUT PARAHETEKS LHAX. HNPIX 
C • AND nXPlX ARE SELECTED 04 T»C Llt^E LENGTH OF THE INKT FILE AND I<ECREE 
C OF ROTATION REOUIHED FCR RECTIFICATICAI. SECTIONS CF THE OUTPUT FILE OF LENGTH 
C LHAX M4D CONTAINING PIXELS ra4PlX TO HXPIX ARE OEN£F;ATED 8Y ONE OR HORE ITER 
C ATIONS OF THIS PROGRAN. 8ECTICR4S Af% CWiCATENATED IN SUDSEOUENT PROCESSING. 

C ACTIONING THE OUTPUT FILE IS DO^E IN RES^SE TO T»« HESSAGEt *1NBI^ ARRAY 
C TOO SHALL." 

C 

C DECLARATICFIS 

REAL NI3).SO> 

LOGICAL*! lNSUF(200000).(»/rB^(3000> 

INTEGER FC8IN(3S>.FCB0UTC^) 

C 

C SET PARAHETEKS 

imiTE(A*99> 

99 FORHAT Cl X. 'ENTER PIXEL/LINE COEFFICIENTS 't > 

READCS.Sfr) R.8 
HR1TE(6.S8> 

88 FORHATtlX. 'ENTER INPUT FIRST PIXEL. NO. PIXELS') 

READ(5.84> IPIXIN.NPIXIN 
HRITEC6.87) 

87 FORHATC IX. 'ENTER HAX. OUTPUT LINES.HIN/HAX PIXEL 'O 

KADCS.84) LHAX.HNPIX.HXP'iX 
84 FORHATCCIA.O) 

c 

C OPEN DATA PILE TO TRANSFORH 

READC8) LENREC. NORECS 

C SETUP WORK ARRAY 

NRECS*niNC 200000/LENREC.NORCCS) 

nXLlNE>M!ECS 

HNLlNE-1 

U1>LENREC*H0DC I . NRECS)*1 

•)2-Jl<»LE»»EC-l 

DO 80 1>I.NRECS 

80 READC8) I1NSUF(J).J*U1.J2) 

C 

C.....OPEN OUTPUT DATA FILE 

mpxoutbhxpi x-ntiF I X* 1 

IFINPXOUT.LE.3000.ANC).NP1X1N.LE.3000>CO to 110 
too STOP ' 7 DATA FILES EXCEED 8CFFER WIDTH.' 
no UR1TEI8) NPXOUT. LHAX 

e.....READ/WRITE LOOP 

DO 210 LINE-1. LHAX 
DO 200 IPX-HNPIX. HXPIX 


9 ^ 


ilil 


PAGE !8 
^ POOR QUALny 


|iNPX-K(t>«(IPX> * R(3)*aiN£» ♦ KO) 
imNPX.i.T.t.OI«.INPX.OT.I#lXlN> THEN 
OUTf^ I l^X'HNP 1 X* 1 )■. FALSE. 

ELSE 

INLlNC*Sm«(IfX> ♦ tlSMILlNC) ♦ SO) 

|f( INL1NE.LT. t. OR. INLINE. CT.NOAECS) THEN 

ouTfuF ( IF i-nm x« i >•. falu. 

ELSE 

C 0«CK LITT FOR SCANLINE IN UORK ARRAY 

IFtlNLlNC .C4.niA.lNE. AND. inline. LE. niLIt4E) 00T0I2O 

C ELSE REAi> SCATA.INE INTO LIST AND WCtRK ARRAY 

IF < INLINE. LT.HNLINE) STOF'INBUF ARRAY TOO SMALL!!' 

00 no l>nXLINE«l. INLINE 
IREC>nODl 1 . NRECS ) •LENRCC 

tlS READ(S.RECbINLINE> tlNDUF(J).J*lREC*l.lREC*NPXIN> 
HXLlNC*lr«.INE 
HNLlNEBMXLlt^-NRECS^l 

120 OUTD>.^(lFX-rwiPlX'»n>lNDUF(LCNREC«KOD(lNLlNC.NRECS>«l!^X) 
ENDIF 
ENDIF 

200 CONTINUE 

210 WRITE(F) (OUTSl^CU>.U*l*»rXOUT> 


FINIS 
200 STOP 
END 




List of Illustrations 


Pigurt l.-^Tha day and night orbital tracks of the HQfl satallita arc inclined 

- ' with respect to esch other, and the area of overlap is a diamond. 

For this example, the scan lines are inclined at an angle of about 
30®. 

Figure 2.~lmagcs of the Powder River Basin area, Wyoming. Figure 2a is the 
transformed night thcratal image, figure 2b is the day theraal 
image, and figure 2c is the original night thermal image. The 
scale of the original images is approximstcly 1:5,000,000. The 
Bighorn Hountalns, Wyoming, arc tc the lower left, the Black Hills, 
South Dakota, are near the bottom center, and the large body of 
water at the top of the scene is Lake Sakakawca, Worth Dakota. 
Significant changes in the appearance of the landscape arc evident 
between the day and night images. The most apparent aspect of the 
affine transformation is a rotation; however, there are also scale 
changes in both directions. The usual convention of light associ- 
ated with high values and dark with low values was employed. 
Extensive cloud cover is present in the upper left h:nd corner of 
the night image and the lower left hand corner of the day image. 


Figure 3 


Figure 4 


.~lMgcs of the Cabeza Prieta arca» aouthweatern United Statea and 
Mexico. Figure 3a ia the day thcraal iaMge, figure '3b ia the night 
thenal tnage, and figure 3c ia a thermal inertia Image constructed 
by applying a model to registered day^thetmal, reflectance, and 
night thermal images. The scale of the original image is approxi- 
mately 1:5,000,000. The Gulf of California is at the bottom of the 
scene and Salton Sea, California is in the upper left corner. 
Clouds are present in the upper right hand corner of the night 
image. 

.'^‘Registration vectors between two sets of control points for the 
powder River Basin area, Wyoming. Figure 4a shows the vectors 

between the day and night images. Figure 4b shows the residual 

vectors after applying the affine transformation. Maximum residual 
is 2 pixels. 
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